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SUMMARY 


A method Is developed for calculating the pressures on aerodynamic 
shapes at" very high supersonic speeds In dense .air -with the ratio of 
specif ic -heats equal to 1. .The method is applicable 'to any body i of 
revolution at zero angle of attack and to any two-dimensional profile. 
The results of the present paper are compared with previous work of . 
von Harman and Epstein .on this, problem and .the differences explained . 

Some aerodynamic characteristics of several shapes are calculated, 
and the lift and drag coefficients are shown to be dependent upon the 
thickness ratio, thickness distribution, and angle of attack. 


INTRODUCTION 


The possibility of constructing airplanes and .missiles capable of 
traveling at extremely high supersonic speeds is of ever Increasing 
importance. Reference 1, for instance, shows that Mach numbers of the 
order of "30 or kO are necessary for obtaining very long ranges with 
rockets . The largest part of such flight paths -will probably be at very 
high altitudes. The aerodynamics of flight at these altitudes has been 
treated by Sanger and Tsien in references 2 and 3* 13a© part of the 

flight in dense air (that is, in a fluid which may be considered a ccm- 
-tinuum) Is also of importance since the surface temperatures and air 
loads may reach very high values. The problem of aerodynamic heating 
has been investigated in ref erence 4. It Is the purpose of the present 
paper to investigate the air loads encountered. 

Busemann has done much of the original work in developing a theory for 
hypersonic flows, however his work on this subject is not generally known 
.in this country. After the present paper had been completed. Dr. Busemann 
made available to the authors a copy of his original work (reference 5 ) • 
The method developed in the present paper Is in substantial agreement 
with that reference. Yon Karmen, in his Yolta Congress address of 1935 
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(reference 6 ), referred to Bueemanr^s theory "but in applying the method 
imposed physically impossible conditions on the flow. Another theory 
has been presented by Epstein (reference 7); however, since he neglected 
the pressure relief due to centrifugal force, his method is applicable 
only to wedges and cones. 
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radius of curvature 

thickness ratio — 

free— stream air velocity 
coordinate - axes 
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0 shock-wave angle 

0 angle "between surface and x-axis 

Subscripts : 

0 free stream; due to form 

a,b ahead of and behind shock wave 

1,2 points of zero pressure coefficient 

c due to centrifugal force 

■i end point of any surface element 

1 lower surface 

L due to lift 


max 


N normal 

s du6 to shock 

T tangential; total due to pressure 

u upper surface 

METHOD OF ANALYSIS 


In order to calculate the pressures on bodies at hypersonic speeds 
it is necessary to determine the position of the nose shock wave. First 
consider two equations from the theory of shock waves. The relation 
between the density ratio across a shock wave and the Mach number normal 
to the shock is 


Pb 7 + 1 



the relation between density ratio and the shock-wave angle 0 and 
deflection angle p is 
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2 - 
If Mjy is very large — — can be neglected and the following- very 

useful approx ima t i on for high Mach number flow is obtained: 


f>b 


r; 


7 + 1 
7-1 


For 7 = 1 

£k->co 

Pa 

and for 7 = l.k - 

_*6 

Pa 

Since 0 is not generally equal to 90° (tan 6 ^ 00 ), f 1 ^ r = 1 

— r ■ 

fc_ 

0 = p - (la) 

and for 7 = lA and slender bodies 

e = 1.20 (lb) 

These relations may be applied only to combinations of free-eatream 
Mach numbers and surface angles which result in very large Mach numbers 
normal to the shock (surface). At the limit M = » aU surface angles 
may be considered. 

The shock equations have been .developed on the assumption of an 
irreversible adiabatic process (no heat transfer except that occurring 
within the shock wave itself). Epstein (reference 7 ) points out that 
the compression at very high Mach numbers results in extremely high 
temperatures which in turn cause large heat _ losses due to radiation and 
conduction. These heat losses limit the temperature rise to a value much 
lower than that determined by the adiabatic law. A better approximation 
to the shock equations at hypersonic speeds may possibly be obtained by 
use of- the assumption that 7=1 in the preceding equations. Equa- 
tions (l) show that the results are not critically affected by the value 
chosen for 7j hence/ for simplicity, the calculations In this paper will 
be based on the relation 7=1. 

Equations (l) indicate - that at very high supersonic Mach numbers 
the shock wave follows the surface. Since the flow Is supersonic the 
air ahead of the shock wave is not affected by the body, and hence the 
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field, of disturbed flow around tie body is limited to a thin layer behind 
the shock which has been called the- hypersonic boundary layer. The 
velocity in this layer may be determined from, the shock equations. 

From the following figure it is apparent that the mass flow (per 
unit area) through the shock ware is 

Pa^%2f - Pb^bjf 



Figure 1.— Velocities near a shock wave. 


Pa 


Since — — =>0 it follows that 
p b 


*N 


u, 


->0 


8-N 


Therefore U- D ^ will be considered negligible relative to U a ^.. 


Equating momentum before and after the shock in a tangential 
direction gives 


(e a ffa]j)Uai - (ptfJbir)^ (3) 

As determined from figure l a equations (l) a (S} a and ( 3 ) the tangential 
Velocity behind the shock (along the surface) is 

U b T = “ U a cos 0 “ cos 0 
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At very high Mach numbers with 7 = 1 the region of disturbed flow 
is confined, to an infinitesimal layer between the surface and the shock 
wave. Since the cross section of this disturbed layer is very small and 
the density is Very large, the acceleration in the direction of flow is 
negligible. Hence the assumption is made that the speed of a given mass 
of the fluid -remains constant along the surface at the value U 0 cos 3. 
This assumption is the same as that made by Busemann and von Harman. A 
sketch of the velocities over a surface is given in figure 2. 



3 sin 3 ^ 


The drag of the body is found from the rate of change of the 
momentum of the air in the stream direction. The air mass per unit time 
entering an elemental area of shock for a body of revolution at zero 
angle of attack is 


2jtrp 0 U 0 dr 

and the change in velocity in the stream direction from the time the 
elemental mass enters the shock at 3 until it leaves the surface being 
considered at 3i is 


U 0 — TJ 0 cos 3 cos 3i 

The drag is then 


D = 



2 *rp 0 ir Q (U 0 - IT C 


cos 3 cos 3 t) dr 
jw v — ' 


00 
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where cos p is a known function of . r, and p^ corresponds to the 
surface angle at the point rq where the flow diverges from the surface. 

The location of rq will.be determined later in the analysis. Equation (10 
differs from, that given' by von Karmen inasmuch as he assumes that the 
flow leaves the surface in the free— stream direction (pq = 0_, rq = r^-y ). 
This assumption would require surface pressures less than absolute zero , 
and is therefore not physically justified. For a two-dimensional body 
the elemental mass is p 0 U 0 dy and the drag for one surface becomes 


D = 




/ P 0 U 0 2 (1 - cos p cos p-j_) dy 

rv 


(5) 


For a cone p remains constant at the same value as Pq. Then the 
drag of a cone at zero angle of attack simplifies to 


r^max 

E = 2ir / P 0 rXJ 0 2 (l - cos 2 Pq) dr 

J 0 


> 


= icc 2 p^U r , 2 sin 2 p 1 
max H o o ^1 


( 6 ) 


or 


Cp = 2 sin^^^ 


(7) 


Equation (6) is identical with Epstein*s result (reference j). Since 
Epstein neglected the effect of the pressure drop due to surface 
curvature , his method would be expected to be in agreement with the 
present work only for straight— side bodies. 

The pressure coefficient behind a shock wave is frequently given as 


2 sin p sin 9 
cos (0 — p) 


(8) 


Equations (l) and (8) combine to give the pressure coefficient behind a- 
shock at hypersonic speeds 


P = 2 sin 2 P 


(9) 
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A comparison of equations ( 7 ) and. ( 9 ) indicates hhair, _ for 7 = 1 , 
two— or three-dimensional "bodies that have no surface curvature in the 
stream direction have the same pressure on the surface as "behind the 
leading shock . "When the flow is curved in the' stream direction, the 
difference in pressure from, the shock wave to the surface equals the 
centrifugal force due to the curvature of the flow. ' The pressure coeffi- 
cient "behind the shock depends on the surface slope, apd the pressure 
relief-due .to curvature depends on the local air mass, velocity, and 
radius of curvature V The maximum curvature that the flow can have for 
a given shock pressure, mass, and velocity occurs when the pressure on 
the airfoil side of the flow drops to zero. If a pressure still lower 
is needed tor turn the flow sufficiently to follow the airfoil or body 
surface, it is physically impossible for the flow to follow the surface 
and it effectively "separates.” For the very high dynamic pressures 
encountered at' hypersonic speeds the pressure coefficient P corresponding 
to zero pressure can be considered as zero. It- follows that the limits 
of integration r^ and y^_ of equations ( 4 ) and (5) correspond to the 

points where the surface pressure coefficients go to zero. Thus, in 
order to determine the limits of integration for the drag integral it is 
necessary to know the local pressure coefficient. 

If in equation (5) the subscripts 1 are replaced by i where 
the subscript i corresponds to the end point of any element, the 
expression will give the drag for the part of the body ahead of the 
point (x^yj.) 


-= P 0 u o 2 f 1 C 1 “ 008 P 008 Pi) 


ay 


In terms of drag coefficient based on unit chord 


r 7± 

Lj_ = 2 / (1 — cos 0 cos ( 3 -^) dy 

‘ u n 


c D; 

The surface pressure coefficient at the point (xj^y^) is then 

ac Dl 


(10) 


Pi = 


ayi 


For ^ x 1 


P-L = 2 sin^ 0 ^ -f 2 sin 


a£i r?i 
ayi L 


cos 0 dy 


(Ha) 
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and for 


z i = *L . t 


*i = 0 


( 1113 ) 


Similarly, for a three-dimensional "body of revolution the elemental 
drag coefficient based on the mirm cross section of the complete body 
is 



r(l — cos p 


cos P^) dr 


and on the surface 


1 dCp^ 

Pl “ as^xi “ 


For x.j_ ^ ' ; ;:1 


P 2 sin P^ 

P., = 2 sin^04 + —= 

i ■ - 1 r. dri 




r i 


r cos p dr 


(12a) 


and for x^ x-j_ 


P ± = 0 (12b) 

■where the fineness ratio F = fe>ng& . Bflga . Jj g. na 

Maximum thickness 

The value of yj_ = y-j_ (or r^ = r^) is determined by equating the 

pressure coefficient to zero. This is the proper value to be used in 
evaluating the total drag integral. (See equations (4) and (5)*) The flow 
separates from the surface at the point (x^y^ ) and is bounded by a 
zero— pressure ' streamline . Any body -within this zfero— pressure region 
contributes no drag. 

In equations (lla) and (12a) the first term is the pressure coeffi- 
cient behind the shock and the second term is the pressure coefficient due 
to centrifugal force. 
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A good approximation to the. surface pressure coefficient for slender— 
profiles may he obtained in a simple manner. Let y = f(x) he the. shape 
of the body or airfoil then 


f(i) 

dx 


and 


dx^. 


f"(x) 


The pressure coefficient behind the shock (equation (9)) becomes 

“ X 



2f ,2 ■ 

1 + f‘2 “ 


If the velocity is not changed appreciably by the shock (true only for 
slender bodies or airfoils}, the pressure drop due to centrifugal 
force is ' ~ 


(Ap) c = 


mTJ o 2 ' 2qy 

b' ” E 


or 



where the radius, of curvature: = : . — 


B 


Ci + r’ 5 ) 3/S 


. f” 


The pressure orr the surface is the shock pressure plus . the change, due 
to centrifugal force: „ : " 

P = P e , + ,P C ■■■■:■■ : - ■ ;■ 


t». . . „ 


2f’ 2 gff 1 

1 + f^ 2 (l + f' 2)3/2 


+ 
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Since f’ 2 is small compared -with 1, the pressure .coefficient on a 
slender airfoil may "be taken as 


P = 2(f.' 2 * ff") 


(13) 


For a three— dimensional tody the air mass over each surface element 


is 


p Q r 

2jtrTJ 0 2 


By similar reasoning the pressure coefficient for slender, three- 
dimensional todies may te approximated as- 

P = 2f‘ 2 + ff" (14) 


The lift of a tody or airfoil results from the downward momentum 
imparted to the air stream. The same considerations used in determining 
the drag from, the horizontal momentum apply to the calculation of lift 
from the -vertical momentum. The lift for the upper surface of a two- 
dimensional airfoil is then of the form • 


L = — p 0 Uo sin Pi 


L 


*i 


cos p dy 


(15) 


where y^_ is the point where the surface pressure "becomes negligible. 


AltaCATION OF THEOKY 


Lift and drag of a flat plate .— Equation ( 9 ) can be used as the 
basis for calculating the "lift of straight— surface airfoils. The lift, 
coefficient of a flat plate at a cana l "1 angle of attack can "be written as 


C L = P Z - P u 

Since the suction pressures are negligible in comparison with the 
pressure rises 

C L = 2<x 2 - 0 = 2a 2 
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and 


Cp — Cj-a — 2a^ 

Lift and drag of -wedge airfoil The lift of a single-wedge airfoil 
is easily obtained by the same method used in the preceding section. For 
angles of attach a less than the semi-nose angle 

C L-=*Z~* U 

C L ~ 2[0 + a) 2 - (3 - a) 2 ] 


and 


— 8j3a = Lta 


dC^ 

da 


- 8p = bt 


(16) 


The drag coefficient of a single-wedge airfoil at an angle of 
attach a which is less than the semi -nose angle is 

Cjjj, = 2 [(e + tt )3 + (e-o) 3 ] 


= 4|3 3 + 12f3a 2 


Ihis can be broken into a form drag C D and a drag due to lift 




pDt, s (17) 


The combination of equations (l 6 ) and (17) giTes the following expression 
for drag due to lift for this airfoil: 

•3 , Crj 2 

Cr. = 4 C T a = - 2 . 

■KL 2 L 3 
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For angles of attack larger than 3, the lift and drag of the -wedge and 
flat plate are the same for the same slopes of the lower surfaces. 

Lift, drag, and pressure on a two— d i mensional -parabolic-a rc airfoil .— 
From equation (5) the drag coefficient for a symmetrical two-dimensional 
airfoil of unit chord may "be expressed as 


C 


D 



(1 - cos 3 cos Pj.) 




For the parabolic airfoil defined by j - itx^l — j!0 where t is the 
thickness ratio, the drag coefficient becomes 


= 2tx 1 (2 - x-|_) 



|/l + t 2 

\]l + t 2 (l — x-^)^ 


(18) 


The value of x^ for which the pressure coefficient is zero is found 
to be 


*1 = 1 - | 1/(1' + t 2 ) 1//3 - 1 


Substitution of x^_ in equation (l8) gives the following equation for 
drag coefficient: 

(1 ; t*) 1/3 _ 1 



or for simplicity of numerical calcuJLations it may be expressed in a 
series as 
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The pressure coefficient on the surfaoe of the airfoil at any* point (xi,yi) 
is found from equations (ll). Then, for x-j_ ^ xq. 


-J: 


/l + t‘ 


\l + t 2 (l - Xi) 2 ] 3 2 j 


and for x* > 


i 


P = 0 


or,, approximately. 


P = t 2 [ 3(1 ~ Xi) 2 - l] + 1 + 6 (1 + Xi) 2 - 15(1 - ij )**! + .. . 

The pressure coefficient found "by the approximate method (equation (13)) 
is, for xq^.x^. 


P = t 2 [ 3(1 ~ Xi) 2 - 1 


The agreement of the pressure coefficients found by these two 
methods is very good as can be seen in figure 3(&) where the pressure 
coefficients on the surface as well as the pressure behind the shock 
are given for a parabolic— er.c airfoil of thickness ratio equal to 0.10. 

The lift coefficient (based on a unit chord) is found from 
equation ( 15 ) to be 


C L = 2 sin 



P dy + 2 sin 



cos P dy 


-*2 


sin p dx + 2 sin P- 


= 2 sin P 2 


o 


o 


sin P dx 
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where (x]_,yj) and are the points of flow separation on the 

upper and lower surfaces, respectively, and J3 is the angle between the 
airfoil surface and stream direction. Let 0 = P — a for the upper 
surface and 0 = ]3 + a for the lower surface where 0 is the angle 
between the airfoil surface and the x— axis . Then for small angles 


sin (0 •— a) 0 — <x 

For the symmetrical parabolic airfoil defined by y 
coefficient based on a unit chord is found to be 

= jt 2 X£>3 — 3tXg 2 (t + co) + 2xg(t + a) 2 J 

-jt 2 :^ - Jtx^Ct - a) + 2x]_(t — a)®j (19) 

The pressure coefficients for the upper and lower surfaces are found 
from, equations (ll). Setting the pressure coefficients found from 
equations (ll) equal to zero permits the determination of the values 
of Xp and Xg 

P u = 3t 2 x^ — 6txp(t — a) + 2(t — a) 2 = 0 


P-£ = 3^ 2 3C2 2 “ 6tX2(t + a) + 2(t + a) 2 = 0 

The lift is then found by substituting these values of x^ and Xg 
into equation ( 19 ). For a <t ~ 


= tx 


(i-l) 


the lift 


2i/3 


°L " ~9 


(t + cl) 3 _ (t — a)3 

t t 


or for t a SSL t 


1 + 1/3 


„ 2/3 

°L " 9 


(t + a)3| 
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or for 



C-^ = t 2 — 3t(t + a) + 2(t + a)^ 


The pressure distribution for a tenr-percent-thick parabolic airfoil 
is shown in figure 3(h) for a = 0.05 radian. If a is larger than the 
semi-nose angle, the effect of the upper surface on the total lift is 
zero since suction pressures are neglected. 

Drag and pressure on a parabolic body of revolution .— The drag of a 
parabolic body of revolution at -zero angle of attack may be determined by 
the same method used for the two-dimensional bodies.- Let the surface be 
described by the equation 




where IP is the fineness ratio of the body. The expression for drag 
coefficient is then • • 


r r 1 

Cp = 1 6 f 2 J (l — cos p cos Pq)r dr 
o 


which reduces to 


= 2x 1 2 (2 - x x ) 2 + 




8f^ 

3 

3 

F 2 


! 

L 


/F 2 + 1 


(1 - x x )2 

52 


- 2 


(F 2 + -l) l / F 2 


F 2 



/I + 


(1 - x^Y 


F 2 


where Xj_ is again the point of zero pressure coefficient,. 
The expression for Cp expanded in a series is 


C D = 


3F 2 


1 + 3(1 - x x ) 2 - 9(1 - *x' + 5(1 - x x ) e 


2hF 


3 + 4(1 - xx) 2 + 18(1 - xx) 4 - 60(1 - xx) 6 + 35(1 - xx) 
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The surface pressure coefficient found from equation (12 ) is, for <^Xj_, 

/ 

2F 2 


P = 2 


1 + 


3xi(2 - x t ) 


<1 ~ 


, i\ 3/2 

1 + w) 


\ 


!+■%(! “ ^ ) £ 


.13/2 


> 


and, for x.j. ^ x^. 


P = 0 


or as a peries, for x^ ^ x^j 


P “ 1 2^(h ^ [' F t 1 - 6(1 - + 5(1 - ^ 

+ ^ [l + 9(! - x ± ) 2 - 45(1 - x ± ) 4 + 35(1 - Xi)^|+ ... > 


(20) 


By use of the approximation given by equation (l4) the pressure coefficient 
may also he expressed, for x^ <L x-j_, as 




5x^ 2 - lOXjL + 4^-^ 


( 21 ) 


The pressure distribution over a parabolic body of revolution of fine- 
ness ratio equal, to 10 as determined by equation (20 ) is plotted in 
figure 4. The value of the pressure distribution as determined by the 
approximate method (equation (2l)) is also plotted in figure 4 together 
■with the pressure behind the shock. 
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. .. CONCLUDING- REMARKS 


A method for calculating the aerodynamic character 1st Ice of two— 
and three-dimensional shapes at hyper sonic speeds in dense air has been 
developed. The results of the present paper are In substantial agreement— 
with some earlier wort of Buaemann but differ from the results of von Kerman 
and Epstein. Ton Kaimian, in his application of Busemann*s method, has 
assumed that the flow leaves the surface parallel to the free stream. 

This assumption would require surface pressures lees than absolute zero. 
Epstein neglected the pressure relief due to centrifugal force and there- 
fore his results apply only to bodies with zero curvature. 

The method of the present paper has been applied to several simple 
shapes and compared with other calculations. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 5, 1948 
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